LOOP SYMMETRIC FUNCTIONS AND FACTORIZING MATRIX 

POLYNOMIALS 



THOMAS LAM 



Abstract. These are notes for my talk at ICCM 2010, Beijing. We survey 
some results, obtained jointly with Pavlo Pylyavskyy, concerning the ring of 
loop symmetric functions. Motivations from networks on surfaces, total posi- 
tivity, crystal graphs, and discrete integrable systems are discussed. 



1. Introduction 

The ring of symmetric functions naturally occurs when studying the problem of 
factorizing a polynomial into linear factors. If one instead considers factorizations 
of matrix polynomials, one obtains the ring LSym of loop symmetric functions. 
Whereas the ring of symmetric functions is the ring of invariants of a natural 
symmetric group action on the polynomial ring, the ring of loop symmetric functions 
is the ring of polynomial invariants of a birational symmetric group action. 

The first half of this article (Sections |[||) is a brief introduction to LSym, and 
certain distinguished elements called loop Schur functions. 

Our work with Pavlo Pylyavskyy on loop symmetric functions began with our 
study of the theory of total positivity of loop groups [LPOS , LP+|. This connection 



is discussed in Section [7[ and the character theory of the infinite symmetric group 
is given as motivation. 

The birational 5 m -action for which LSym is the invariants arises as the birational 
iJ-matrix of certain affine geometric crystals [BK, Et, KNC| . We explain this 



connection at the level of combinatorial crystals in Section |8|, and show how loop 
Schur functions can be applied to the study of the energy function of affine crystals 



LPlOa . The energy function (for the affine crystals in this article) is equivalent to 



Lascoux-Schiitzenberger's cocharge statistic on tableaux, and it is the connection 
to the latter that we describe. 

This birational symmetric group action also leads to discrete integrable systems. 
We describe Takahashi-Satsuma's box-ball system [TS] in Section^, and explain the 



relation, due to Hatayama, Hikami, Inoue, Kuniba, Takagi and Tokihiro [HHIKTT 



In this context LSym can be viewed as integrals of motions of these dynamical 
systems. 

In Sec tion p|, we explain the connection between LSym and networks on the 
cylinder [ LPlOcf . On the one hand, this theory is a powerful technique for estab- 
lishing properties of the birational symmetric group action. On the other hand, 
these topological networks lead to a broad collection of generalizations. 

The birational 5 m -action of this paper also arises in work of Noumi and Yamada 



NoYa on discrete Painleve systems, though we shall not discuss this connection. 
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My results reported upon in this article are all joint work with Pavlo Pylyavskyy, 
and I am grateful to him for our long collaboration. 

2. Symmetric polynomials 

2.1. Via factorizing polynomials. Let pit) = (1 + x\t){\ + x 2 t) •••(! + x m t) 
be a polynomial with roots — x^ 1 , —x^ 1 , ■ ■ ■ , — x^ 1 . The elementary symmet- 
ric polynomials e 1 {x 1 ,x 2 , ■ ■ .,x m ), e 2 (x 1 ,x 2 , ■ ■ .,x m ), . . ., e m (x 1 ,x 2 , ■ ■ -,x m ) are 
defined by the equation 

p(t) = l + e 1 t 1 +e 2 t 2 + --- + e m t m . 

Thus 

eifzi, . . .,x m ) =x± + x 2 H hi™ 

e 2 (xi, . . . , x m ) = xix 2 + xix 3 H 1- x\x m + x 2 x 3 H h x m ~ix m 

and so on. 

2.2. Via an action of the symmetric group. Let the symmetric group S m act 
via algebra automorphisms onZ[i 1 ,i 2 ,..., x m ] by the formula 

(1) W ■ f(xi,X 2 , ■ ■ .,X m ) = f{x w (x),X w {2),. . .,X w ( m )). 

The following "fundamental theorem of symmetric polynomials" is often attrib- 
uted to Newton. 



Theorem 2.1. The elementary symmetric polynomials ei,e 2 ,...,e m are alge- 

■ • ) Xrn\ 



braically independent generators of the ring of invariants 1[x\, x 2 , . . . , x r , ] " 



We shall denote the ring of Theorem 2.1 by Sym m , and call it the ring of 
symmetric polynomials in m variables. The inverse limit of the Sym m , is the 
ring of symmetric functions in infinitely many variables, denoted Sym. 

2.3. Via Galois Theory. The symmetric group S m also acts on the field Q(xi,x 2 , ■ ■ ■ , x m ) 
of rational functions in m indeterminates via the same formula (|l|) . In this context 
one may think of S m as the Galois group of the Galois extension Q(xi, x 2 , • ■ ■ , x m )/Q(ei, e 2 , . 
which has degree ml. 

This point of view has the weakness that one has to work with fields, rather than 
rings. 



3. SCHUR POLYNOMIALS 



We refer the reader to [EC2| for more concerning the material of this section. 
The ring Sym m has a Z-basis consisting of the Schur polynomials sx(xi, x 2 , . . . , x m ), 
where A = Ai > X 2 > • ■ ■ > A m > is a partition with at most m parts. 

3.1. As a generating function of tableaux. If p C A as Young diagrams, then 
we can define the skew Schur polynomial 



sx/fj,(xi,x 2 , . . .,x m ) = 2J 



x wt(T) 
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where the sum is over all skew semistandard Young tableaux (SSYT) with shape 
A/yU, filled with the letters 1,2, ... ,m and wt(T) denotes the weight of the tableau 
T. For example, 









1 


1 


3 


5 


■5 


T = 


1 


2 


2 


2 


5 


G 




4 


4 


6 













is a SSYT with shape (8, 7, 3)/(3, 1) and x wt < T ) = xfx^xjxlxj. When fi = (0) 
then s\/(o) — s x is a Schur polynomial. 

3.2. As a Jacobi-Trudi determinant. If A is a partition, then A' = (A' l7 . . . , A^) 
denotes the conjugate partition, where A • is the number of boxes in the i-th column 
of A. 

Theorem 3.1 (Jacobi-Trudi formula). The skew Schur polynomial S\/^ has the 
following determinantal expression in terms of elementary symmetric polynomials: 

s\/n = det(e A /_ i+j _^.)™. =1 . 
In particular, the skew Schur polynomials Sx/^xi, . . . , x m ) lie in Sym m . 

3.3. As a ratio of alternants. Let a = (a\, a 2 , . . . , a m ) be a sequence of non- 
negative integers. We define the alternant 

a «= E (-!r™(z Q ) 

wes m 

where (— l) w denotes the sign of w. The alternant can also be expressed as the 
determinant a a = det(x" 3 Define the staircase partition 5 m := (m,m — 

1,...,0). 

Theorem 3.2. For A with less than or equal to m parts, 

, \ _ n A+g m _i 

S X {-^1 1 • • • 7 %m ) 



4. Loop symmetric polynomials 

Fix an integer n > 1, which will usually be suppressed in our notation. We shall 
define a ring LSym m of loop symmetric polynomials in m sets of variables. 

When n = 1, the ring LSym m reduces to Sym m . 

4.1. Via solving matrix equations. Suppose one is given a n x n matrix P(t) 
with coefficients which are polynomials in t. Is it possible to factorize P(t) into a 
product of "linear" matrix factors? In the following we will formulate precisely one 
version of this problem. 

Define the n x n whirl matrix 

1 fli ••■ \ 
1 a 2 ■■■ 
1 ••• 



M(a 1 ,a 2 ,...,a n ) 



I) (j 

\ a n t 



1 a„_i 
1 J 
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For n — 1, we define Ad (a) = 1 + at. We think of M(ai, d2, ...,a„) as a polynomial 
(in i)-valued matrix, which in particular is linear. Every matrix entry of M has 
degree less than or equal to 1. 

Now let P(t) be a n x ji matrix with polynomial coefficients. Can one factorize 

(2) P{t) = M{x<p , . . . , 4"» ) • • • M(x$ , . . . , xt ] ) 

(k) 

for suitable "roots" x\ ? There is a natural condition to enforce: each matrix 
M(x^~\ . . . , il 1 ') becomes an upper triangular matrix with l's on the diagonal 
when t — 0. So to obtain a solution, we must at least assume that P(0) is an upper 
triangular matrix with l's on the diagonal (uni- upper triangular). This is analogous 
to the fact that the constant term of pit) = (1 + xit){l + x 2 t) ■■■{! + x m t) must be 
1. 



Theorem 4.1 ([LPlOb, Theorem 4.1]). Suppose P(t) is anxn matrix with entries 
in C[t], such that P(0) is uni-upper triangular. If P(t) is generic, then it has a 
factorization of the form (||) . 

Given a factorized matrix polynomial P(t) = M{x^\ . . . jX^) ■ ■ ■ M(xm\ ■ ■ ■ ,Xm^), 
we define the loop elementary symmetric functions {e$ | i = 1, . . . , n and k = 
1, . . . , m} by the matrix equation 

/ l + + + + eg>_ a t + ■ • ■ J» 1+I ^_ 1 t + ... \ 

P(t) = 



4- i) t+ei" + - 2 v + ... e ("- 1) t+ e «T 3 1) * 2 + e<"- 1) +e, < r + " 1 1) *+- 

V 4 n) * + 4+i* 2 + ■ • • 4" ) *+4+2* 2 + 4-i* + 4"-i* 2 + • • ■ i+4 n) t+--- 

For example, let us take n — 2 and m = 3. One has 



r (i) T (2) (l). 




so that 



Jl) _ r (l) , r (l) , r (l) p (2) _ _(2) , (2) , (2) 

„(1) _ r (l) T (2) , T (1) T (2) , _(1)_(2) (2) _ (2) (1) (2) (1) (2) (1) 

(_^2 — 1 "-"2 I «*•< 2 3 I «// J-- g ^2 — 1 2 3 I ^2 3 

_ _(1)_(2)_(1) _(2) _ (2) (1) (2) 

3 — 1 2 3 3 — 1 2 3 

Note that is a homogeneous polynomial of degree i. An explicit formula for 

(r) . 

(r), \ ( r ) ( r+1 ) (r+fc-l) 

e k \ X li ^2) ■ ■ • , ^mj — / 4 x i 1 x i 2 ' ' ' x i k 

l<ii <%2< - ■ ■ <ik<m 
(r) (r) 

By convention, el = for A; < 0, and eg — 1. We call the upper index the color. 

When all n colors are identified, that is x[ s ^ — x\ s ' for all i and s, s' € Z/nZ, these 
polynomials specialize to the usual elementary symmetric polynomials. We denote 
the subring of Zfa;^ , x^ 2 ' , . . . , x^ , x^ , • ■ ■ , im'] generated by the loop elementary 
symmetric polynomials e^ by LSym m , the ring of loop symmetric polynomials 
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in m sets of variables. Our loop elementary symmetric polynomials (with some 
change in labeling of variables) appeared previously in work of Yamada 



Remark 4.2. In \LP0^ we distinguished two different rings of loop symmetric 



functions, one called the curl ring, and the other the whirl ring. The two are 

(r) (—r) 

related by negation of the upper index x\ >-> x\ , where as usual the upper 
index is taken modulo n. 

4.2. Via the birational i?-matrix. The whirl matrices M{x) do not in general 
commute with each other for n > 1, but they do satisfy a commutation relation. 
For two n-tuples (o^ 1 ', . . . ,x^ n ') and (i/ 1 -*, . . . ,y^ n ') of variables, we denote 



i+n— 1 j i-\-n— 1 

*(*,»)= e n y (k} n 

j—i k—i+1 k—j+1 



where here and elsewhere, upper indices are always taken modulo n. For example, 
for n = 3, one has 

K 1 (x,y) = x^x^+y^x^+y^y (3 K 
Now define a rational map 

a : Q0r« . . . , x (n \y (1 \. . . , y™) -> . . . , xW.jW, . . . , yW) 

by 

(3) fl ( g (<)) =y (<+9 g!±feg) and S (y«) = ^-l)^^ . 

Ki( x ,y) K i{ x 7y) 

If we are given to sets of variables {(x^, . . . , xj n ^) 1 < i < m}, then we let 
s/c denote the rational map acting on ®(Xj), fixing the sets of variables Xi for 
i 7^ k, k + 1, and then applying s to the variables x^, Xk+i- 

Theorem 4.3. 

(1) The rational maps Sj generate a birational action of S m on ). 

(2) Let w G S m act on Q(x^) as in (1). Then 

M{ Xl )M{x 2 ) ■ ■ -M(x m ) = M{w{x l ))M{w{x 2 )) ■ ■ ■ M(w(x m )). 

Note that if n = 1, then this birational action of S m is just the usual action 
of S m on Q(xi, . . . , x m ) of Section |^. The birational action of S m , and the proof 
of Theo rem Q(l) was established independently in a number of different contexts 
El], g |NoY4 |LP08| [LPlOcf . Theorem pj|(2) is phrased in terms of "M-matrices" 



in the context of affine geometric crystals 



KNO|]. 



The following result will be established in [ LP+ . Its proof is significantly harder 



than the fundamental theorem of symmetric functions. 

Theorem 4.4 (Fundamental theorem of loop symmetric functions). The ring of 
loop symmetric polynomials is exactly the polynomial invariants of the birational 
S m -action: 

LSym m = Q(xf ) ) s <»nZ[xf ) ]. 

Furthermore, the loop elementary symmetric polynomials are algebraically indepen- 
dent generators o/LSym m . 
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We note that the definition of LSym m is compatible with restriction of variables: 
setting x m = gives LSym,,,^ from LSym m . Thus it is possible to define a ring 
LSym, in infinitely many variables, as the inverse limit of LSym m . We call LSym 
the ring of loop symmetric functions. 



Example 4.5. Let n — 3 and m — 2, and set x — x\ and y — X2- Then 

s(xW +y (1) ) 

; (2)^ + 2 .(3)^3 



y 



x (2) x (3) + y(2) x (3) +y (2) y (3) 



,(*) 



,(!) 



— -"1 tI/j — c x . 

4.3. Galois groups of matrix polynomials. It can be deduced from the pre- 
ceding discussion that Q( X \^)/Q( X 'f ^) Sm is a Galois extension of degree m!, with 
Galois group S m . (Indeed, this is significantly easier than Theorem 4.4, see [LPlOb, 
Theorem 4.1].) Thus one can think of S m as the Galois group of the matrix polyno- 
mial P(t), where the coefficients of the polynomial are the "variables" e^p . Perhaps 
it would be interesting to ask: what can be said about the "Galois groups" of matrix 
polynomials with coefficients taking values in Z or Q? 



5. Loop skew Schur functions 

5.1. As a generating function of tableaux. A square s = in the i-th 

row and j-th column has content c(s) = i — j. We caution that our notion of 
content is the negative of the usual one. Let p/v be a skew shape. Recall that a 
scmistandard Young tableaux T with shape p/v is a filling of each square s G p/v 
with an integer T(s) G Z>o so that the rows are weakly-increasing, and columns 
are increasing. For r G Z/rtZ, the r- weight x wt< ( T ) of a tableaux T is given by 



x 



.wtC»(T) 



n 



s&p/v *T(„) 



(c(a)+r) 



We shall draw our shapes and tableaux in English notation: 



o o o o o 
o o o 
o o o 





1 


1 


1 


3 


1 


2 


2 


3 


4 




3 


3 


4 









For n = 3 the 0-weight of the above tableau is (x\ ') (x\ ) x\ x 
We define the loop (skew) Schur function by 



(1)^2^(1)^3^(2)^(2)^,(2) (3) (3)/ (3h 2 
i x l ^2 \. x 4 ) ■ 



.wtW(T) 



where the summation is over all semistandard Young tableaux of (skew) shape X/ p. 
Example 5.1. Let n = 2. Then 



(1)/ ^ _ (1) (2)„(2) 

(1) (2) (2) 
*^ 1 1 

corresponding to the tableaux 



, r (l) r (2) T (2) , 



UU 2 ) (2) 



1 


1 


2 





1 


2 


2 





1 


2 


3 





(1) (2) (2) 
X i X 2 



(1)^(2)^(2) 



1 


3 


2 





T (1) T (2) T (2) 
(1)^(2)^(2) 
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1 


1 


3 





1 


3 


3 





where the upper left corner has content 0, and so gives a 1-weight with color 1. 
Setting x^ — xf^ — Xi gives the usual Schur polynomial S2,i(xi, x 2, X3). 

Perhaps surprisingly, the loop Schur functions do not form a basis for LSym. By 
Theorem 4.4, monomials in the loop elementary symmetric functions will form a 



basis for LSym, but finding the analogue of monomial symmetric functions is not 
straightforward. The following problems seem central. 

Problem 5.2. Find a Schur-like basis for LSym. 

Problem 5.3. Find a monomial-like basis for LSym. 

Presumably any "loop monomial symmetric function" would be nonnegative in 
terms of monomials of the variables, but be minimally so in some sense. 



5.2. As a determinant. We have the following analogue of the Jacobi-Trudi for- 
mula. 



Theorem 5.4 ( LP08, Theorem 7.6]). The loop skew Schur function has the fol- 
lowing determinantal expression in terms of loop elementary symmetric functions: 



(r) 



det(e£:^; 



Thus 



£ LSym. 



The proof of Theorem 5.4 is essentially the same as that of Theorem B.l 



Remark 5.5. The Jacobi-Trudi formula puts loop Schur functions into the context 
of Macdonald's ninth variation of Schur functions. The latter have a similar Jacobi- 



Trudi formula NNS_Y], for which the upper index is not cyclic. 



5.3. As a ratio of alternants. Suppose we are given annXm array {x+ | 1 < 
i < m and j € Z/nZ} of variables. Let r 6 Z/nZ. Given a decreasing sequence 
a.i > a 2 > • • • > a n > of nonnegative integers, we define the loop alternant 

(r) 

act &S the m x m determinant 



det i t -,i (V( r M r_1 ) ■ ■ • T ( r - a i+ 1 ) 



)rn 



where t a j, denotes the transposition exchanging a and b, acting via the birational 



action of S m . The following result [LP+ generalizes Theorem B.2 



Theorem 5.6. For A with less than or equal to m parts, the loop Schur function 



lr-l) 



can be expressed as 



l X+5 



In an earlier version of this paper, this result was stated as a conjecture. We thank Greg 
Anderson for the idea of the proof. 
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Example 5.7. Let A = (2, 1), n = 3 and m = 2, denoting the variables by ir = xi 
and y = X2- Then 

(i) _ j f y (1 V 3 V 2) X W X W X <2) 



«a = det( ^ x(3), ; 3/Kl 

So 

_ t, (i y 2) !/ (3) x (3) (i (1) x (2) + j/ 1 M 2 > +2/ (1 V 2) ) -lOiWilVlfiOid + w (2) x (3) + y (3) y (3)) 



3.(1)3.(2)3.(3) - y(l)y(2)y(3) 



(3) 



5.4. LSym as a Hopf algebra. The ring of symmetric functions has a Hopf al- 
gebra structure, where comultiplication is given by A(ej) = Y^j=o e 3 ® e i-h an< ^ 
the antipode is given by S(ei) = {—l) l hi, where hi is the complete homogeneous 
symmetric function. We now briefly describe the Hopf structure of LSym. 
The comultiplication of LSym is given by 

A(ef>) = X>f 0e^> 
3=0 

where by convention = 1 for every k. Note that LSym is not co-commutative. 
The antipode is given by 

Clearly LSym is graded, with deg(e[ k ^) = i. 

Proposition 5.8. The above structures define a graded Hopf algebra LSym. 



Let MSym, the Hopf algebra of matrix symmetric functions be the graded 
dual Hopf algebra to LSym. By general theory, one knows that MSyniQ = MSym(g>z 
Q is generated (as an algebra) by its primitive subspace prim. In fact MSyniQ is 
the universal enveloping algebra of the Lie algebra prim. 

Proposition 5.9. The primitive subspace prim of MSyniQ is the the Lie-algebra 
of infinite Z x Z indexed matrices, which are n-periodic, strictly upper-triangular, 
and which contains finitely many non-zero diagonals. The grading on prim is such 
that the i-th diagonal above the main diagonal has degree i. 

5.5. Loop powersum symmetric functions, and a loop Murnaghan-Nakayama 
rule. For k > 1, define the loop powersum symmetric polynomials 

Pk{x\, X2, • • • , x m ) = (x\ ^x\ ^ ■ ■ ■ x\ ^) + (x 2 ^ x 2 ' • • • x 2 ^) + • • • {x[^ x[^ ■ ■ ■ x^) 

so that pk is homogeneous of degree kn. The Pk{xi, x 2 , ■ ■ .) are primitive elements 
of LSym, though even the fact that they lie in LSym is not completely obvious. 
Recall that a ribbon is a connected skew shape not containing any 2x2 square. 
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Theorem 5.10 ( jLP+| ). We have 

Pk 



4 r) = E(-D 



ht( M /A) (r) 



where the summation is over all ribbons /i/X of size kn, and ht([i/\) denotes the 
number of rows in jJi/X, minus 1. 



Thus the usual Murnaghan-Nakayama rule | EC2 | holds for the loop powersum 
and loop Schur functions. 

6. Networks on cylinders and the birational S^-action 
6.1. Local description of S^-action. The fact that the transformation (H) satis- 



fies the braid relation (Theorem 4.3(1)) is not obvious. Here we sketch an approach 
based on a topological model of networks developed in [ LPlOc ], 

Let us first arrange our variables xf ^ on a network embedded into a cylinder. The 
variables are located on the vertices of this network, that consists of n horizontal 
wires connecting the boundaries of the cylinder, and m loops which go around the 
cylinder. The case n = m = 2 is illustrated in Figure |l|. The birational action 
s : Q(x\^) — > Q(xp' ) ) will be expressed in terms of certain local transformations of 
the network. 





A 






A 










! i i x'M 








V 


V 


i \ i 


V 


V 



Figure 1. The birational action (x, y) (x',y'). 



We need the following two local transformations (more are discussed in [LPTOc]), 
the Yang-Baxter relation (Figure^) and the crossing creation/removal rule 
(Figure |3|). The parameters in this Yang-Baxter rule naturally occur in Lusztig's 
study of total positivity in reductive groups |Lus| . 



t x + z 



Figure 2. Yang-Baxter move with transformation of vertex 
weights shown. 

To calculate the birational action s, we first use the crossing addition rule to 
create a pair of crossings with weights p and —p. Using the Yang-Baxter move, 
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Figure 3. Adding or removing a crossing. 



we may push one of these crossings through to the horizontal wires, until it loops 
around the cylinder. In effect, this has swapped the two vertical loops. For a unique 
choice of the weight p, the weight that comes out on the other side after passing 
through all the horizontal wires is also p; and so one can use the crossing removal 
rule the pair of crossings with weights p and —p. 




Figure 4. Swapping two vertical loops by creating a pair of cross- 
ings, pushing one of them around the cylinder, and removing these 
crossings. 



Theorem 6.1 ([LPlOc]). The above local procedure calculates the birational S n 
action of Section (.!t. 



Using Theorem 5.1, one can prove many properties of the birational S' m -action, 
including Theorem |4.3| (1), using "local" computations. 

6.2. Boundary and cycle measurements. Take a network N m such as the one in 
Figure [j], but with m vertical loops and n horizontal wires (as usual n is suppressed 
from the notation) . Then the loop elementary symmetric functions can be obtained 
by taking certain generating functions of paths in this network. 

Every intersection of these networks consists of two wires intersecting transver- 
sally. Suppose p is a path. Then at each intersection it can either go straight, or 
turn. A highway path is a path which always turns when the crossing wire crosses 
your current wire from the left. The weight of a highway path is the product of 
the vertex weights over all vertices where the path goes straight. See Figure |[ 

Lemma 6.2. Loop elementary symmetric functions are exactly weight generating 
functions of highway paths in N m with fixed starting and ending boundary vertices, 
and fixed homology class. 

The generating functions of the above lemma are called boundary measure- 
ments. 

Example 6.3. Let us take the network in Figure [I]. Then there is only one highway 
path which starts at the upper left vertex and ends at the upper right vertex, with 
weight e% (x, y) = x^y^ 2 \ This is the highway path which just goes straight. 
There are two highway paths from the upper left vertex to the lower right vertex, 
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Figure 5. The three ways a highway path can go through a vertex 
with weight x; a fragment of a highway path, contributing the 
factor xz to the weight of the whole path. 



with weights x^ and j/ 1 - 1 respectively, summing to e^ 1 "*. 
paths loop around the back of the cylinder. 



Both of these highway 



In Example 3.3, we did not have to fix the homology class of the path. But we 
would have to do so for larger m. 

One can also dehne cycle measurements. In the case of the networks N m , we 
should consider underway cycles, which are cycles in N m where the condition at 
each intersection has left and right swapped, compared to highway paths. Underway 
cycle measurements in N m give the loop powersum symmetric functions. (There is 
a technical multiplicity which is involved in the definition of cycle measurements, 
which we omit here.) 

Example 6.4. In Figure |l|, there are two underway cycles which loop around the 
cylinder exactly once. These have weights x^'x^ 2 ' and y^y^ respectively, which 
sum to p\{x,y). If we consider underway cycles which loop around the cylinder 
twice, we will get two cycles, with weights (x^x^) 2 and (j/ 1 ^ 2 ') 2 , summing to 

P2(x,y). 

The fact that boundary and cycle measurements are preserved by the birational 
SVn-action is not a coincidence. In [LPlOc , we establish 



Theorem 6.5 (jLPTOcJ Theorem 3.2] 
preserved by local transformations. 



Boundary and cycle measurements are 



Here "local transformations" consist of a list of moves, which includes those in 
Figures || and [| 

Remark 6.6. The topological interpretation of LSym gives rise to many gener- 
alizations. In [ LP10c| , we also prove that one obtains an action of the symmetric 
group when the horizontal wires do not point in the same direction (see Remark 3.3 
for an application of this). However, any (valid) network on an orientable surface 
would give rise to a ring of invariants, where the symmetric group is replaced by a 
"monodromy group" of networks we define in [LP 10c . 



7. Total positivity 

7.1. Total positivity implies factorization. A polynomial p(t) with complex 
coefficients always has a factorization into a product of linear factors. The analogous 
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statement is not true for the matrix factorizations 

„(!) 



P(t),¥(^,...,,H-W,...,i"») 

we consider, but there is an elegant condition which always guarantees factorization 
over the reals. 

Given a matrix polynomial P(t), we associate to it an infinite periodic matrix 
X = Xp, as follows: 

/ 



1 + At + 9t 2 
4 + 6i 2 
4 



+ t 







2-2* 
7 + t-t 2 
-2 
1 



+ t 2 





-1 















4 -2 
1 


9 

6 -1 










1 2 
4 7 


4 -2 
1 


9 
6 -1 










1 2 
4 7 


4 -2 
1 















P(i) 



X P 



In other words, Xp is the block- Toeplitz matrix, where each block is equal to the 
coefficient of t l in P(t), for some i. 

Suppose P(t) is a matrix polynomial with real coefficients. We now define P(t) 
to be totally nonnegative if every minor of Xp is nonnegative. 



Theorem 7.1 ( ILPOSH ). Suppose P(t) is a real matrix polynomial, with P(0) uni- 
upper triangular. Then P(t) is totally nonnegative if and only if there is a factor- 
ization of P(t) into a product of whirls with nonnegative real parameters. 

For n = 1, we recover a well-known criterion for a polynomial to have nonnegative 
real roots. 

A much more difficult problem is to classify the real matrices with formal power 



series coefficients which are totally nonnegative. This is studied in [DMS, LP08|. 
This problem can be reformulated in terms of LSym. Let f7>o denote the set of 
matrices P(t) with real formal power series coefficients such that P(0) is uni-upper 
triangular, and such that P(t) is totally nonnegative. 

Proposition 7.2. The set U>q is in bijection with the set 



{4> : LSym <g) 

of algebra homomorphisms of LSym i 
all loop skew Schur functions. 



>R\cb(s[ r l)>0} 

to the reals taking nonnegative values on 



When n = 1, the set U>o is the set of totally positive functions: formal power 
series p(t) = 1 + a%t + a^t 2 + ■ ■ ■ such that the infinite Toeplitz matrix 



/ 1 



a i 
1 



«2 

ai 
f 



«3 

ai 
f 



is totally nonnegative. 
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7.2. Relation to infinite symmetric group. The classical Frobenius character- 
istic map identifies the ring Sym of symmetric functions with the direct sum of the 
character rings of the symmetric groups {S n \ n > 1}. It is perhaps not surpris- 
ing that the set of totally positive functions has an interpretation in terms of the 
representation theory of the symmetric groups. 

The infinite symmetric group Soo is the union U n >iS n of all the finite symmetric 
groups, where S n -\ is identified naturally with a subgroup of S n . A character of 
Soo is a complex-valued function x '■ 5oo — > C which satisfies four conditions: 

(1) central - x is constant on conjugacy classes 

(2) positive-definite - the matrix ix{9^ 9j))\ j=i 1S Hermitian and nonnegative- 
definite for any g u g 2 , . . . ,g r 

(3) extremal - x is n °t a sum of two linearly independent functions satisfying 
(1) and (2) 

(4) normalized - %(1) = 1. 

For a finite group, such a function would exactly be an irreducible character, nor- 
malized by dividing by the dimension of the character. 

Let us call a totally positive function p(t) = 1 + ait + a^t 2 + ■ ■ ■ normalized if 



a% = 1. The following result is due to Thoma [Th] 



Theorem 7.3. The set of normalized totally positive functions is in bijection with 
the set of characters of the infinite symmetric group. 



Via this theorem, the homomorphisms of Proposition 7.2 have the following 
interpretation. Given a character x '■ &oo C we may restrict to S„ to get a 
(reducible, fractional) character x\s n of S n . The value <fi(s\) of the corresponding 
homomorphism <f> : Sym — > K is the coefficient of the irreducible character \ X m 

X\s n - 

If a totally positive function p(t) is a polynomial, then it factors into linear factors 
of the form (1 + at). In the case that p(t) is not a polynomial, the Edrei- Thoma 



theorem [Edr ( Th] gives a complete description of p(t) as a meromorphic function. 
In this context, Vershik and Kerov JVKj give a beautiful interpretation of the poles 
and zeroes of p(t) : thinking of a character of as a point-wise limit 

X = lim X XM /x Xin) (l) 

n— >oo 

of normalized irreducible characters of S n , the poles and zeroes control the asymp- 
totic growth of the lengths of rows and columns in A*™). 



While the following problem was one of our main motivations in [LPOS], little 
progress has been made. 

Problem 7.4. Find an interpretation of totally nonnegative matrix formal power 
series (U>o for n > 1) in terms of representation theory of Soo, or of some other 
asymptotic representation theory. 



As Proposition |7.2| suggests, this may require expressing LSym as the character 
ring of a series of groups, in the way that Sym is the character ring of the series of 
symmetric groups. 

8. Crystals, energy, and charge 



In Section 4.2, we discussed a birational action of S m on a n x to array of 



variables. This action of S m turns out to be the i?-matrix of certain products of 
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Kirillov-Reshetikhin crystals. We first give an informal introduction to Kashiwara's 
crystal graphs; focusing in particular on the relation with tableaux. 

8.1. AfRne i?-matrix and jeu de taquin. For a combinatorial introduction to 



this subject, we recommend Shimozono's "Crystals for dummies" Sh]. Crystal 
graphs were invented by Kashiwara [Kas| as combinatorial models for representa- 
tions of quantum groups. For our purpose, a crystal graph (of affine type A) is a 
finite graph with edges labeled by one of 0, 1, 2, . . . , n — 1. Here is an example of 
an affine crystal graph: 



1 


1 


1 


1 


2 


l 


2 


2 








For a positive integer s, the set of semistandard Young tableaux of row shape 
(s), with entries in 1, 2, . . . , n, form an affine crystal, which is a special case of a 
Kirillov-Reshetikhin crystal. 

If £>i, B 2 are Kirillov-Reshetikhin crystals, the combinatorial i?-matrix is the 
unique isomorphism Rbx.b 2 '■ B\ ® B% ~ ► B 2 ® B\ of affine crystals. It is known 
that the combinatorial R- matrices generate an action of S m on B\ ® • • • ® B m . 

The i?-matrix has a convenient combinatorial interpretation ph| in terms of 
semistandard tableaux and the jeu de taquin algorithm [EC2|. Let b\ ® b 2 be an 
element of B\ <g> B 2 . Then Rb u b 2 (bi <8> b 2 ) = c\ ® c 2 S B 2 ® B\ where c\,c 2 are the 
unique pair of row shaped tableaux which jeu de taquin to the same tableau that 
b\ and b 2 jeu de taquin to, as follows: 





1 


1 


3 


3 


3 


1 


2 


3 







1 


1 


3 


1 


2 


3 


3 


3 





since both jdt to 



1 


1 


1 


3 


3 


3 


2 


3 





The action of t he combin atorial i?-matrix can be described in terms of coordi- 
HHIKTTfl. If x = 



! X 



— 



Jn) 



,.T 



('<- 



nates, as follows 
x the n-tuple (x^ 

subtraction-free rational function / is obtained by replacing 
+, and / by — . 



) is an n-tuple, we denote by 
Recall that the tropicalization trop(/) of a 

by min, and x by 



(r) 

IJ2 be the 
. , n . Then 



Theorem 8.1. Let RB 1 ,B 2 (bi ® 62) = c\ ® C2 and let x^, x% , yi~\ 
number of boxes filled with r-s in b\, 62, Ci, C2 respectively, for r = 1, . 

trop(s(xi,x 2 )) = (3/1,2/2)- 

Example 8.2. In the example above, we have x\ = (1, 1, 1), x 2 = (3,2,0), y\ — 
(1, 1, 3) and y 2 — (1, 2, 0). To calculate y\ , one takes 



Vi = trop(x 2 ' 



(3) K3- 
K 2 ' 



min(2, 4, 5) - min(2, 1, 3) = + 2 - 1 = 1. 
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8.2. The energy function and charge. The energy function Db ■ B — >• Z of an 

affine crystal plays an important role from two different points of view: (1) from the 
representation theory point of view the energy encodes information corresponding to 
the affine weight of elements of the crystal, and (2) in the connection between affinc 
crystals and vertex models, elements of a tensor product of crystals are thought of 
as paths, and the energy function has the physical interpretation as the energy of 
a path. See [ |KKMMNN ] . 

Instead of giving the usual axiomatic or recursive definitions of the energy func- 
tion, we will instead describe it in relation to another well known statistic on words: 
cocharge. Suppose u is a word in the letters 1,2, ... ,n. Then the weight wt(u) of 
u is the composition a = {a\, . . . , a n ) where a, is the number of i's in u. Suppose 
the weight of u is a partition, namely, ai > ot2 > • • • > ot n . The cocharge cc(u) 
is obtained recursively as follows. Underline the rightmost 1 in u. Given an un- 
derlined i — 1, underline the rightmost i to its left, if it exists; otherwise, underline 
the rightmost i in u, if it exists; and otherwise stop. To each underlined letter 
we assign an index: the 1 is assigned index 0, and the index of the underlined i 
is equal to the index of the underlined i — 1 if the i occurs to the left of the i, 
otherwise the index of the i is one more than the index of the i — 1. Finally, we 
define cc(u) — indices + cc(u') where the word v! is obtained from u by erasing 
all the underlined letters. We define the cocharge of the empty word to be 0. 

Example 8.3. The cocharge of u = 3222311111233 is calculated as follows (indices 
are indicated as subscripts): 

3022203111110233 

220311110233! 

2 31 11023! 

3illo2i 

lo 

Thus cc{u) = (0 + + 0) + (0 + + 1) + (0 + + 1) + (0 + 1 + 1) + = 4. 



Let T be a semistandard tableaux with partition weight [1%, . . . , jj, m . Let 
(T) be the number of i's in the j-th row of T. We define a crystal element 



b{T) - h 



1 b m £ B b 



B„ 



tableau with i's. These b(T) are exactly the highest weight vectors of the affine 
crystal. Recall that the reading word r(T) of a tableau is obtained by reading the 
rows from left to right, starting from the bottom row. The following result is due 



B^ m by setting hi to be the one-row 



to Nakayashiki and Yamada [NoYa]. 



Theorem 8.4. Let T be a semistandard tableaux with partition weight. Then the 
cocharge cc(r(T)) of the reading word is equal to the energy Db(1)(T)) of the corre- 
sponding crystal element. 

Example 8.5. Let 





1 


1 


1 


1 


1 


2 


3 


3 


T = 


2 


2 


2 


3 












3 
















Then r(T) is the word in Example 


3.3. Thus the 


energ 



b(T) 



1 


1 


1 


1 


8) 


1 


2 


2 


2 


8> 


1 


1 


2 


3 
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is 4. 



It is known that the energy function of a tensor product i3 Ml (g> 



)Bu_ of 



Kirillov-Reshetikhin crystals commutes with the affine i?-matrix. As explained in 



Section 4.2, the ring of loop symmetric functions is exactly the ring of polynomial 
invariants of the birational i?-matrix. Nevertheless, the following result is somewhat 
surprising: 



Theorem 8.6. Set x 

<g> ® • • ■ (g 
(o) I 



a,- J . The energy function Db of the affine crystal 
B fim is equal to the tropicalization of the loop Schur function 



■ j x m ) . 



There seems to be no a priori reason that the energy function Db should be the 
tropicalization of a polynomial, much less that of a loop Schur function. 



Example 8.7. For Example 3.5 we would have X\ 
X3 = (1,1,2). According to Theorem [ 



(5,0,0), x 2 = (0,1,3), and 



the energy is 



min(xj 



(•'!) 



X 



X 



(3) 



' ] ~^ *^1 

corresponding to the tableaux 



(2) 
1 

(2) 



X 



X 



(1) 
1 

(1) 



X 



X 



(3) 

1 

(3) 



4 1} 



,.(3) 



T (3) T 0) 
x 3 1 x l 



,(3) 



T (l) 

x 2 



„(3) 



r (l) 

x 2 



1 


1 


1 


1 


2 


2 





1 


1 


2 


3 


3 


3 





1 


2 


2 


2 


2 


3 





The minimum is achieved (uniquely) on the tableau 



1 


1 


2 


3 


2 


3 





(3) , J 2 ) , J 1 ) , J 3 ) . (1) . ^(3) 



which gives x[ 
with the previous calculation. 



= + + + 2 + + 2 = 4, agreeing 



9. Box-ball systems 

Loop symmetric functions can also be thought of as integrals of motion of certain 
discrete dynamical systems, called box-ball systems. Here the birational action of 
S n is the time evolution of the system. The connection between the i?-matrix and 



the box-ball system was established in [HHIKTT 



9.1. Takahashi-Satusuma box-ball system [TSJ. We begin with an (infinite) 
configuration of boxes arranged on a line. Each box can either contain one ball or 
no balls. We shall always assume that there are finitely many balls. Time evolution 
proceeds as follows. Take the leftmost ball and place it in the leftmost empty box 
to the right of that ball. Repeatedly do this for the leftmost ball that has not yet 
been moved, until all balls are moved. This completes the time evolution. 

Figure ^ shows such a configuration, and its time evolution. 

While the definition is simple, the box-ball system has the remarkable property 
that it exhibits solitonic behavior. An isolated consecutive string of k balls is a 
soliton which travels with speed k. Figure ^| shows the collision of a speed 3 
solition with a speed 1 soliton. Note how the faster soliton overtakes the slower 
one, but both solitons are preserved. 
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9.2. Carrier description, and the birational i?-matrix. The time evolution 
of the box-ball system can also be described in terms of a carrier. Imagine that a 
carrier (with infinite capacity) begins from the left going from box to box. At each 
box, the carrier: 

(1) Picks up a ball, if the box has a ball in it. 

(2) Drops a ball, if the box is empty and the carrier is carrying a ball. 

It is easy to see that this carrier description gives the same outcome as the time 
evolution described previously. Let us rewrite the interaction of the carrier with a 
site in terms of coordinates. A site with a extra capacity and b balls will be denoted 
(a. b); thus the empty box is denoted (1, 0), and the box with a ball in it is denoted 
(0, 1). The carrier will be denoted (oo, d) if it has d balls in it. So when the carrier 
begins from the left, it is in the state (oo, 0). Then the interaction of the carrier 
and the box can be written as 

{(a, b), (oo, d)} i — > {(a — min(a, d) + 6, min(a, d)), (oo, d — min(a, d) + b)}. 

The reader is invited to verify that this transformation is (after some relabeling) 
the tropicalization of the birational ^-action for n = 2: 



FIGURE 6. Time evolution in the box-ball system. 




It follows that Theorem 4.4 can be interpreted as a statement concerning the 



integrals of motion of the box-ball system. 
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Proposition 9.1. The polynomial integrals of motion of the birational analogue of 
the box-ball system is exactly LSym. 

If we imagine that the system has m boxes, then the invariants "total number of 
balls" and "total amount of free space" are equal to the tropicalizations of e^' and 

(2) 

e m . Because of the curious way the box-ball coordinates are shifted with respect 
to the LSym-variables, the invariant e± can be given the following interpretation. 
Its tropicalization is equal to 1 if every odd site is empty and every even site is 
occupied, and equal to otherwise. Each time evolution, the notions of odd and 
even swap (pushing the carrier through shifts the parity). It would be interesting 
to give simil ar interpretations for a general loop (skew) Schur function. 

See TTSfl for a description of invariants of the box-ball system in terms of the 



Robinson-Schensted correspondence. 

Remark 9.2. There is a vast literature generalizing box-ball systems to situations 
with many different kinds of balls (corresponding to our n > 2), or to situations 



where boxes have higher capacity. See [HHIKTT| for some directions. 



Remark 9.3. Together with Pylyavskyy and Sakamoto, we give in [LPS| a gen- 
eralization of the box-ball system which arises from considering mixed directional 
horizontal wire on a cylinder, as in Section [j[ This new solitonic system involves, 
boxes, balls, and baskets. (Baskets can contain balls, but they are allowed to move 
as well.) 
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